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Abstract 

We study two-magnon spectra in resonant inelastic x-ray scattering (RIXS) from the Heisenberg 
antiferromanets by extending the formula of Nomura and Igarashi (Phys. Rev. B 71, 035110 
(2005)). The core- hole potential in the intermediate state of RIXS gives rise to a change in the ex- 
change coupling between 3d electrons, leading to an effective interaction between the core hole and 
spins of 3d electrons. We derive a formula suitable to calculate the two-magnon RIXS intensities, 
replacing the bare core-hole potential responsible to charge excitations by this effective interaction 
creating two magnons in our previous formula. It consists of two factors, one of which determines 
the incident-photon-energy dependence and another is the two-magnon correlation function. We 
evaluate the former factor for La2CuC>4 in terms of the density of states of the 4p states obtained 
by the band calculation. We also calculate the two-magnon correlation function as a function of 
energy loss u> and momentum transfer q of the Heisenberg model on a square lattice, by summing 
up the ladder diagrams after transforming the magnon-magnon interaction into a separable form. 
The calculated spectra form a broad peak around uj = 3 J for S = 1/2 on the magnetic Brillouin 
zone boundary and vanish at q = (0,0) and (it, 7t). Such momentum dependence of the RIXS 
spectra could provide an excellent opportunity to study the dynamics in the Heisenberg model. 

PACS numbers: 78.70.Ck, 72.10.Di, 78.20.Bh, 74.72.Dn 
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I. INTRODUCTION 



Resonant inelastic x-ray scattering (RIXS) has recently attracted much interest, since it 
provides valuable information of charge excitations in solids.— >2»M&&iI Unlike optical mea- 
surement, it can probe directly the momentum dependence of the excitations. The i^-edge 
resonance in transition-metal compounds is particularly useful, because the corresponding 
x-ray wavelengths are an order of lattice spacing. In this situation, the Is core electron is 
prompted to an empty 4p state by absorbing photon, then charge excitations are created in 
order to screen the core-hole potential, and finally the photo-excited 4p electron is recom- 
bined with the core hole by emitting photon. In the end, charge excitations are left with 
energy and momentum transferred from photons. 

For analyzing such spectra, Nomura and Igarashi (NI) have developed a formalism&*^ 
by adapting the resonant Raman theory of Nozieres and Abrahams.— According to the 
formula NI have derived, the RIXS intensity is expressed in terms of the density-density 
correlation function in the equilibrium system under the Born approximation to the core- 
hole potential. Describing the electronic states on the multiband tight-binding model within 
the Hartree-Fock approximation, and taking account of the electron correlation within the 
random phase approximation, Nomura and Igarashi have analyzed the RIXS spectra of 
undoped cuprates.-&i£ The calculated spectra have reproduced well the experimental ones 
as a function of energy loss and the dependence on momentum. The use of the Born 
approximation to the core-hole potential has been examined by evaluating higher-order 
corrections, and have partly been justified in spite of a strong core-hole potential.— In 
addition, the RIXS spectra in NiO have recently been analyzed by the same method, in an 
excellent agreement with the experiment.— Therefore, the NI formula seems quite useful to 
analyze the RIXS spectra. Note that most theoretical studies on the momentum dependence 
of the RIXS spectra have been based on the numerical diagonalization method for small 
clusters with replacing the 4p band by a single level.— ^ By these numerical methods, it 
seems practically impossible to analyze the RIXS spectra in three dimensional systems. 

Quite recently, Hill et al. have reported that the RIXS intensity has been observed around 
the energy loss 300 ~ 600 meV in L^CuC^ with improving the instrumental resolution.— 
One scenario for the origin of the spectra is that the intensity arises from two-magnon 
excitations. In this paper, we examine this possibility by developing the NI formalism.-"^ 
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As pointed out by van den Brink,— the presence of the core-hole potential in the intermediate 
state modifies the exchange process, giving rise to the change in the exchange coupling 
between the spins of 3d electrons at the core-hole site and those at the neighboring sites. 
This leads to an effective interaction which creates two magnons by the core hole. Replacing 
the bare core-hole potential responsible to charge excitations by the effective interaction 
in the NI formula, we immediately obtain the formula of the two-magnon RIXS spectra. 
It consists of two factors, one of which gives the incident-photon-energy dependence and 
another is the two-magnon correlation function. The former factor involves the density 
of states (DOS) of the 4p states, and is almost independent of energy loss. We evaluate 
it for La2Cu04, using the 4p DOS given by the local density approximation (LDA). A 
large enhancement is predicted at the K edge as a function of incident-photon energy; the 
enhancement is much larger for the polarization along the c axis than for the polarization 
in the ab plane. 

The two-magnon correlation function determines the dependence on the energy loss to 
and the momentum transfer q. We calculate it within the first order of 1/5* (S is the 
magnitude of spin) at zero temperature on the antiferromagnetic Heisenberg model in a 
square lattice. It is known that the 1/5* expansion works well with taking account of the 
quantum fluctuation. The linear-spin-wave theory is made up of a leading-order in the 1/S 
expansion,—*^ and the magnon-magnon interaction as well as the Oguchi correction to the 
magnon energy arise in the first order of 1/S.— ^ Various physical quantities such as the 
spin-wave dispersion, the sublattice magnetization, the perpendicular susceptibility, and the 
spin-stiffness constant, have been calculated up to the second order in 1/5* in the square 
latticed 22 ' 23 ' 24 - 25 In this paper, the magnon-magnon interaction is taken into account by 
summing up the ladder diagrams with transforming the interaction into a separable form. 

It is found that the spectral shape as a function of energy loss is strongly modified by the 
interaction. Such an effect is known in light scattering, although the momentum transfer is 
limited to zero.— i2L22i22i22i2L RIXS could detect the momentum dependence of the spectra. 
Our results show that the spectra form a broad peak around uo = 3J (J is the exchange 
constant) for S = 1/2 at the magnetic Brillouin zone boundary and vanish at q = (0, 0) and 
(7r,7r). Note that the momentum dependence of two-magnon excitations has been studied 
in the context of the phonon-assisted photon absorption spectra,— and has recently been 
discussed in light scattering.^ The latter case may be relevant to the non-resonant inelastic 
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x-ray scattering experiments, but the mechanism causing two-magnon excitations is quite 
different from RIXS without core hole. This may make the RIXS spectra quite different 
from the light scattering spectra because of the different matrix elements. 

The present paper is organized as follows. In Sec. [TTJ we formulate the RIXS spectra 
for the Heisenberg antiferromagnets by adapting the NI formalism. In Sec. 1111} the RIXS 
formula is expressed with magnon operators. The two-magnon correlation function is calcu- 
lated within the first order in the 1/S expansion. In Sec. IIV| numerical results are presented 
for the RIXS spectra in a two-dimensional Heisenberg model for L^CuO^ Section |V] is de- 
voted to the concluding remarks. The spectra (at q = 0) in light scattering are summarized 
in Appendix for the square lattice. 



II. TWO-MAGNON PROCESS IN RIXS 



At the transition-metal K edges, the ls-core electron is excited to the 4p band in a dipole 
transition by absorbing photon. This process may be described by 

H x = wJ2 -4= E < a) Pl^c m e^ + H.c, (2.1) 

qct V H jrja 

where represents the 77-th component (rj = x, y, z) of two kinds of polarization vectors 
(a = 1,2) of photon. Operator c qQ , stands for the annihilation operator of the photon 
with momentum q and polarization a. Since the Is state is so localized that the Is — > Ap 
transition matrix element is well approximated as a constant w. Annihilation operators Pj va 
and Sj a are for the 4p v and Is electrons with spin a, respectively, at the transition-metal 
site j. The Hamiltonians for the core electron and for the 4p electrons are given by 

His = e u ^2 s )a s i^ ( 2 - 2 ) 
H ip = ^eykjp^p^. (2.3) 

In the intermediate state, the attractive core-hole potential works on the 3d electrons, which 
may be described by 

imaa' 
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Annihilation operator di a is for the 3c? state with spin a at site %. Finally, for describing the 
low energy behavior, we assume a single band Hubbard model for 3d electrons, 

H = t J2 + H - c -) + u E d \Ai d ^- ( 2 - 5 ) 

<i,j> i 

Here U may be 4 — 8 eV, smaller than V. This Hubbard model may be mapped from a more 
precise "d-p" model for cuprates. 

At the half-filling, a spin singlet pair has the energy 2t 2 /U lower than a spin triplet pair 
has. Therefore, in the low energy sector, the system may be described by the Heisenberg 
model with the exchange coupling constant J = 4t 2 /U. At the core-hole site, this exchange 
process may be influenced by the core-hole potential, as shown in Fig. (TJ which leads to the 
change of the exchange coupling. This has been pointed out by van den Brink.— The energy 
difference between the spin triplet and singlet of two electrons, one at the core hole site and 
the other at a nearest neighbor site, is estimated as 



t 2 



1 1 

+ 



(2.6) 



u-v u + v_ 

The first and second terms arise from the process that two electrons are on the core-hole 
site as shown in Fig. QJb) and that two electrons are on the nearest neighbor sites to the 
core-hole sites, respectively. Taking the difference of the energy from the value without the 
core hole, we obtain an effective interaction between the core hole at site n and spins of 3d 
electrons, 

Hhs-3d = S ncrSl lcT S n • S n+ S, (2-7) 

a,5 

with 

U 2 V 2 

Jc = Uu 2 ~^T 2 ' (2 ' 8) 

where n + 8 represents a nearest neighbor site to site n. 

The RIXS intensity is simply given by replacing Hi s ^ d with H^f_ 3d in the NI formula.- 
The corresponding diagram in the Keldysh scheme is shown in Fig. [2, where the effective 
interaction H^_ 3d is represented by the double line. The incident photon has a momentum 
qj and an energy Ui, a polarization e^° l \ and the scattered photon has a momentum and 
an energy Uf, a polarization e'*A The momentum and the energy transferred from photon 
are given by q = q^ — q^ and u = Ui — Uf. They are written simply as q = (q, u). Similarly, 
we introduce the notations g« = (q f , and qf = (q^Uf). 



(a) (b) (c) (d) 



FIG. 1: A schematic representation of the two-magnon process in RIXS, corresponding to the 
first term of Eq. (|2.6p ; (a) the Is electron is excited to the 4p band by absorbing x-ray, (b,c) an 
exchange process takes place under the influence of the core-hole potential, and (d) the 4p electron 
is recombined with the ls-core hole. 

The upper triangle represents the product of Green's functions of the 4p electron and the 
core hole on the outward time leg, which gives a factor, exp[i(e4 p (p) — eu — iTi s — Ui)t], with 
Tis being a lifetime broadening width of the Is core hole. The lower triangle represents the 
product of Green's functions on the backward time leg, which gives a factor exp[— i(el p (p) — 
6i s + iTis — Ui)(t' — u')]. Note that extra time dependent factors, e luls on the outward time 
leg and e~ luJS ' on the backward time leg, arise from the interaction. Integrating the time 
factor combined to the above product of Green's functions, with respect to s and t in the 
region of t < s < 0, — oo < t < 0, we obtain 

L n B {ui]to) = J c J dt— ^exp[i(e^ p (p) - e Xs - iT u - uJi)t] J dse luJS 

= _ Jc / ^ (£)d£ . (2.9) 

J (ui + eis + iT Xs -e){uji-uj + e Xs + iT ls - e) 

Here the summation over the momentum is replaced by the integration over the energy 
associated with the 4p DOS pL(e) projected onto the r] (= x, y, z) symmetry. The integration 
with respect to s' and t' on the backward time leg gives the term complex conjugate to 
Eq. (12.91) . The integration with respect to u' gives the energy conservation factor, which 
guarantees that to in Eq. (I2.9p is the energy loss, to = tOi — uj. See Ref. Q for the detailed 
derivation of this function. Since the magnon energy is usually an order of 0.5 eV, which is 
much smaller than the energy scale of the 4p DOS, we can safely put u = in Eq. (12.91) . 
The bubble with shaded vertexes in Fig. [2] denotes the two-magnon correlation function 
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FIG. 2: Diagrams for the RIXS intensity in the Keldysh scheme. The double lines represent the 
effective interaction between the core hole and spins of 3d electrons. The solid lines with 4p and Is 
denote the Green's functions of the 4p electron and of the Is core hole, respectively. The bubble 
with shaded vertexes stands for the two-magnon correlation function, which connects the outward 
time leg on the top half and the backward time leg on the bottom half. 



with 



Y + -(q,uj) 



(Mt( s ')M q ( S ))e^ s '- s )d( S '- S ), 



n,5 



(2.10) 



(2.11) 



where (• • • ) indicates the average over the ground state. Combining all the above factors, 
we finally obtain the expression of the RIXS intensity as 

2 



Iwl 4 N 

W(q h af, q f , a f ) = - — Y + (q, u) 



uj e 



(2.12) 
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III. THE 1/S EXPANSION 

We carry out systematically the l/S'-expansion by introducing the Holstein-Primakoff 
transformation to the spin operators.— Assuming two sublattices in the AF ground state, 
we express spin operators by boson operators as 



Sf 




(3.1) 


sf 


= (Sr)t = V2Sh{S)a u 


(3.2) 


si 


= -S + bfy, 


(3.3) 




= (Sry = V2Sb]f 3 (S), 


(3.4) 



where aj and bj are boson annihilation operators, and 

with m = a\di and bjbj. Indices i and j refer to sites on the "up" and "down" sublattices, 
respectively. 



A. Heisenberg Hamiltonian 

First we apply the 1/S" expansion to the Heisenberg Hamiltonian described by 

H = J S * ' S i' ( 3 - 6 ) 

<i,3> 

where indicates the sum taken over nearest neighbor pairs. Substituting Eqs. fl3.1|) - fl3.4p 
into Eq. (13.61) we expand the Hamiltonian in powers of 1/S as 

H=~JS 2 Nz + H + H 1 + --- , (3.7) 

where N and z are the number of lattice sites and that of nearest neighbor sites, respectively. 
The leading term Hq is expressed as 

H = JSJ2 ( a l a i + h ) h 3 + a ih + a l b ])- ( 3 - 8 ) 

<i,j> 

Rewriting the boson operators in the momentum space as 

/ 2 \ 1/2 

ai = ( /vJ ^2a k exp(ik ■ r 4 ), (3.9) 
/ 2 \ 1/2 

h i = [n) E 6kex P^ k - r ^' ( 3 - 10 ) 



we diagonalize Hq by introducing the Bogoliubov transformation, 

a l = 4«k + fn k p_ k , 6_ k = m k al + 4/?-k, 

where 



2e k J 



1/2 



m k 



1-ek 
L 2e k J 



1/2 



-x k 4 



with 



£k 



1 - 7k> 7k 



k-<5 



(3.11) 
(3.12) 
(3.13) 



Here <5 is the nearest neighbor vectors. Momentum k is defined in the first magnetic Brillouin 
zone (MBZ). 

As shown in Ref. |2l|, after the Bogoliubov transformation, the Hamiltonian becomes 



(3.14) 



JSz 
2S 



A^e k (4ak + At#k) 



1234 

,-: < 



+ (2a t 1 /?_ 2 a3«4 J Bi 4 2 ) 34 + 2/?l 4 /3_ 1 /3_ 2 « 3J Bg 34 + aJa^Is/^aL + H -c) 



with 



A 



N 



CkJ- 



(3.15) 



(3.16) 



The first term in Eq. (13.151) arises from setting the products of four boson operators into 
normal product forms with respect to magnon operators, known as the Oguchi correction.— 
For the square lattice, A = 0.1579. The second term represents the scattering of magnons. 
Momenta ki, k 2 , kj, ■ ■ • are abbreviated as 1, 2, 3, ■ ■ • . The Kronecker delta <5g(1 + 2 — 3 — 4) 
indicates the conservation of momenta within a reciprocal lattice vector G. Explicit expres 
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sions for B^'s in a symmetric parameterization are given by Eqs. (2.16)-(2.20) in Ref. 

(3) 

Here we only write down the explicit expression of B 12U , which will become necessary in 
the next section, 



B 



(3) 
1234 



= 72-4 + 7l-3XlX 2 X 3 X 4 + 71-4X1X2 + 72-3^3X4 
+ 73X2X3X4 + 74-2-lXl + 73-2-IX1X3X4). 



(3.17) 



B. Two-magnon operator 

Inserting Eqs. (13.11) ~ (13. 4p into Eq. (12. lip , we expand M q in terms of boson operators 

as 



Note that the momentum transfer q is defined in the first BZ, which is the double of the 
first MBZ. When q is outside the first MBZ, it can be brought back to the first MBZ by 
a reciprocal vector Go, that is, q = q + Go with q being inside the first MBZ. In this 
situation, e iq ' rj = — e iq °' rj in the second term of Eq. (13.181) . Noting this fact and substituting 
Eq. (13. lip into Eq. (13.181) . we obtain the expression of M q as 



N(q,k) = fe{(l±7 qo )^ qo+k] m k + sgn(7 G )(±l + 7 qo )m [qo+k] 4 

+ [7[q +k] ± sgn(7 G )7 k ]m [qo+k] m k + [ 7k ± sgn(7 G )7tq +k]]^[q +k]4}- (3-20) 

Sign ± corresponds to the case that q is inside or outside the first MBZ. Symbol [q + k] 
stands for the q + k reduced in the first MBZ by a reciprocal lattice vector G, that is, 
[q + k] = q + k — G, and sgn(7 G ) denotes the sign of 7 G . This result will be used to 
calculate Y + ~(q, uo) in the next section. Note that iV(q, k) = for q = and q = (jr, tt) in 
the square lattice, indicating that no RIXS signal would be generated. 

C. Two-magnon correlation function 

Defining the two-magnon Green's function as 




(3.18) 




(3.19) 



with 




(3.21) 



we rewrite the two-magnon correlation function Y + (q, u) as 



F + -(q, u) = |p E N & k ^ N ^ k ') (-2)ImF(q , w; k, k'). 



(3.22) 



k k' 
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Here T is the time-ordering operator. The two-magnon Green's function is expanded in 
terms of the one-magnon Green functions, 

G aa {k,t) = -z(T[e* k (*)4(0)]), (3.23) 

G W (k,t) = -z(T[/3 k (04(0)]>. (3.24) 

The unperturbed ones corresponding to H are given by 

G®(k,u) = &$(k,u) = [u-e k + i7 } ]- 1 , r^0 + , (3.25) 

in the energy unit of JSz. Hereafter we use this energy unit. 

In the lowest order, the two-magnon Green's function is simply given by 

F(q , U] k, k') = F (q , w; k)5 k y, (3.26) 

with 



F (q ,^;k) = % J G^([q + klco + k )G^(-k,-k ] 



dko 
27 



(3.27) 



w - e [qo+k] - e k + it] 

Inserting this relation into Eq. (13.221) . we obtain the correlation function in the lowest order 

as 

y(°) + -(q, u) = Yl N (& V 26 ^ - %o+k] " Ck). (3.28) 

k 

In the first order in 1/5, the magnon energy e k is changed into (1 + ^)e k due to the 
Oguchi correction. Therefore, F (q, u;;k) is modified by replacing q qo +k] and 6[ k ] in Eq. 
(I3.27P with (1 + ^|)e[q +k] and (1 + ^)e k , respectively. Let -F (q, u;k) be the function 
after this modification. Within the same order, we have to take account of the magnon- 
magnon interaction. Only the terms of a factor -Bi 2 34 i n Eq. (13.151) is relevant to the present 
calculation. Consider the ladder approximation shown in Fig. [3] for the two-magnon Green's 
function. Regarding the dependence on k and k' in F(q ,u;k,k') as a matrix F(q ,uj) 
with y x y dimensions, we notice that the sum of the ladder diagrams is equivalent to the 
expression, 

^(q ^) _1 ] kk , = ^o(q ,^; k)- 1 ^' + ^^[ qo+ kl4A qo+ k]4'5g ( ) )+k , ]iki[qo+kLk ,. (3.29) 
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(a) 

a 

[qo+k] 



a 

[qo+k] V 



[qo+k'] 



-k' 



a 
[qo+k] 



a , 
[qo+k'] 



(b) 03 



X/ " 4 r< 3 > iitt X/ 
~2SN 1234 1 2 3 4 



CCl 



IB- 



FIG. 3: (a) Ladder diagrams for F(q , u; k, k'). Solid lines represent the single-magnon Green's 
function with the magnon energy including the Oguchi correction, (b) Magnon-magnon interaction 
in units of JSz. 

This is nothing but the eigenvalue equation for two-magnon excitations, indicating that 
the ladder approximation together with the Oguchi correction to the single-magnon energy 
constitute the first-order correction in the 1/S expansion. 

Equation (13.291) is not useful for the actual calculation of the two-magnon Green's func- 



tion, because the matrix with 



dimensions has to be inverted in order to get F(q , u) 



We sum up exactly the ladder diagrams, transforming the interaction into a separable form 
with several channels, 



4 



2SN 



mil Vn(4, 1). 



^3.30) 



m,n=l 



Here N c is the channel number. The indices 2 and 3 specify the incoming magnons while 
1 and 4 specify the outgoing magnons (Fig. H2(b)). Explicit forms of v n (k,k') and T mn are 
given for the square lattice in the next section. Thereby we obtain the T-matrix LT, 

II(q , u : k, k') = A J2 ^(k, [q + k])rf n (q , u)v n (k', [q + k']), (3.31) 



where 



with 



[I-fi^qo,^]-^ 



2 x -> 

[R(q , u)] mn = j^2^v m % [q + k])F (q , u; k)v n (k, [q + k]). 



(3.32) 



(3.33) 
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In Eq. (13.321) . the unit matrix 1 and R(q ,u) are in the N c x N c dimensions. We calculate 
the two-magnon Green's function from the T-matrix, 

F(q , u; k, k') = F (q , w; k) [5^ + n(q , u : k, k')F (q , w; k')] . (3.34) 

Inserting this equation into Eq. (13.221) . we obtain y +_ (q, u;), which gives the RIXS intensity 
as a function of momentum and energy transferred from photon. 

IV. CALCULATED RESULTS 

We apply the formulas in the preceding sections to La 2 Cu0 4 , which seems to be a typical 
two-dimensional Heisenberg antiferromagnet. The RIXS intensity in the two-magnon region 
relative to that in the charge excitation region is roughly estimated as (J c /V) 2 . For La 2 Cu04, 
we have ( JjVf < 0.01, because J c ~ 0.2 - 0.5 eV and V ~ 5 - 10 eV. 

A. Incident-photon-energy dependence 

We carry out the band calculation for La2Cu04 within the LDA. Figure H^a) shows the 4p 
DOS projected onto the p z and p x (p y ) symmetries, which may correspond to the absorption 
coefficients for the photon polarization along the c axis and in the ab plane, respectively. 

Using the same DOS, we calculate the factor determining the incident-photon-energy 
dependence, | ^ e?, LP^uj^uj = 0)e^| 2 , for the incident and scattered photons having the 
same polarization. Figure IH(b) shows the calculated results for the polarization along the 
c-axis and in the afe-plane. A strong resonant enhancement is predicted as a function of 
incident photon energy. The enhancement for polarization along the c axis is much stronger 
than that in the ab plane. This is quite different from light scattering where the polarization 
of light is restricted in the ab plane (see Appendix).— 

B. Momentum dependence as a function of energy loss 

We consider the Heisenberg model on a square lattice. We put 7k = (cosk x + cosk y )/2 
with k x , k y in units of l/(lattice constant). For the magnon-magnon interaction, we have 
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10 20 
co[eV] 



30 



FIG. 4: (a) The Ap DOS projected onto the p z and p x (p y ) symmetries, calculated by the LDA 
in La2Cu04. The origin of the energy is the bottom of the Ap band, (b) The RIXS intensity 
proportional to | e ^ ^(ujf, u> = 0)e^\ 2 , with e^ along the c axis and in the ab plane. The 
origin of the energy corresponds to the photon energy of exciting the Is core electron into the 
bottom of the Ap band. 



the channel number N c = 13, and f is given by 



r = 



2SN 
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TABLE I: Definition of the coefficients v n (k, k') 



n 


^(k,k') 


1 


i44/x k 


2 


1 44' cos ^ 


3 


1 44' cos ^3/ 


4 


|44' x k' cos(A; a ; — 


5 


^44' x k' cos(k y — k' y 


6 


f44' x k£k' cos k' x 


7 

8 


|44' x k 2; k' cos 

1 44' sin fcx 


9 


f 44' s i n 


10 


|44 ,a V sin^ - fc^ 


11 


l44/x k / sin(A; y - fc^ 


12 


f44' x k^k' s i n ^a: 


13 


f44' x k^k' s i n 



The v n (k, k')'s are defined in Table [B In the actual calculation, we divide the first MBZ 
into 512 x 512 meshes, and carry out the sum over k in Eq. (13.331) to evaluate .R(q ,u;). 

Figure \5\ shows the calculated RIXS intensity F + ~(q, uj) scaled by (Sz) 2 as a function of 
energy loss u for several typical values of q. The lowest order value, F^ + ~(q, u)/(Sz) 2 , is 
independent of S, corresponding to S = oo. The effect of the magnon-magnon interaction 
becomes weaker with increasing S, and the spectra approach to the lowest-order values. As 
already pointed out, no RIXS intensity comes out at the T point and at q = (n,7r). When 
q is close to the T point (panel (a)), a sharp peak is found for u < JSz, which is slightly 
modified in the presence of the magnon-magnon interaction. When q is at the boundary of 
the first MBZ (panels (b) and (c)), a sharp peak found in the lowest order approximation is 
smeared out to be a broad peak due to the magnon-magnon interaction. The center of the 
shape remains nearly the same after taking account of the interaction. In contrast to these 
cases, when q is outside the first MBZ (panel (d)), a sharp peak found in the lowest order 
approximation is changed into a broad peak with its center considerably shifted to the lower 
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FIG. 5: The RIXS intensity Y^ (q, u) scaled by (Sz) 2 as a function of lo/JSz for typical values 
of q. The bold and thin solid lines stand for Y H (q, lo)/{Sz) 2 with 5* = \ and 1, respectively. The 
broken lines represent Y^ ^ (q, ui)/(Sz) 2 , which correspond to S — > oo. 



energy region due to the magnon-magnon interaction. 

Figure [6] shows the RIXS spectra as a function of energy loss with changing momenta 
along several symmetry lines (S = 1/2). Along the zone boundary of the first MBZ (panel 
(b)), the spectra have large widths around uj = 3 J. The spectra obtained here seem to 
be different from the results of van den Brink,— who made the moment analysis. We find 
no long-lived virtual bound state of two magnons, which has been predicted on the phonon 
assisted light absorption spectra.— Note that the matrix elements of creating two magnons in 
RIXS are different from those in light scattering.— Since the spectral shape depends strongly 
on the matrix elements, the direct comparison of the two spectra may be less useful. The 
present formalism recovers the results of light scattering (at q = 0) by Canali and Girvin,— 
as shown in Appendix. 

The exchange coupling constant J in La2Cu04 is estimated as J = 135 meV by comparing 
the spin-wave velocity calculated in the first order of 1/S with the experiment.— Therefore, 
the broad peaks are located around uj = 400 meV for q at the boundary of the first MBZ, 
Unfortunately, no experimental data are available on such a low energy region in I^CuO^. 
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FIG. 6: The RIXS intensity F + ~(q, uj)/(Sz) 2 along the symmetry lines for q within the first order 
of 1/S = 1/2). 

V. CONCLUDING REMARKS 

We have formulated the two-magnon spectra of RIXS in antiferromagnets by developing 
the formalism of Nomura and Igarashi. The Is core hole potential causes a change in the 
exchange coupling between 3d electrons, resulting in the two-magnon excitation. This is 
analogous to the conventional RIXS process that the charge excitation is created due to the 
screening of the core hole potential. The intensity of two-magnon RIXS is estimated to be 
lass than 0.01 of the intensity coming from the charge excitation. 

In the present formalism, the factor describing the incident-photon-energy dependence is 
separated from the factor describing the dependence on the momentum and the energy trans- 
ferred from photon. We have calculated the former factor using the 4p DOS for L^CuC^ 
within the LDA. We have predicted a strong enhancement of the intensity at the K edge for 
the polarization along the c axis. The latter factor is given by the two-magnon correlation 
function. We have calculated the correlation function up to the first order of 1/S in the 
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square lattice, systematically applying the 1/S expansion. We have exactly summed up 
the ladder diagrams by transforming the magnon-magnon interaction into a separable form 
with 13 channels. The spectral shape as a function of energy loss is strongly modified by 
the magnon-magnon interaction. On the boundary of the first MBZ, for example, the sharp 
peaks found in the lowest-order approximation have considerably been broadened. We hope 
the spectra obtained in this paper would be compared with the experimental data in future. 

No AF-long-range order could exist at finite temperatures in purely two-dimensional 
Heisenberg models. In the absence of the AF order, however, it is known from the non- 
linear u model analysis^ that the spin-spin correlation length is rather long up to T ~ J/ks 
The spin-wave-like excitations could exist in such a situation.— Therefore, the present 
analysis of the RIXS spectra at zero temperature may have relevance to the spectra at finite 
temperatures. The analysis of temperature effects is left in future study. 
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APPENDIX A: TWO-MAGNON LIGHT SCATTERING 

We summarize the two-magnon excitation of light scattering to compare the result by 
the present formalism with previous studies in two dimensions. 

Since the wavelength of light is much longer than the lattice spacing, the intensity is 
independent of the directions of the incident and scattered lights. The interaction of light 
with spins is described by^ 

H R = £ [A(et ■ S)(e f ■ 6) + B(e t ■ e f ) + Ce*e* f ] (S j ■ S j+S ), (Al) 
jS 

where 8 is a unit vector in the direction joining the nearest neighbor pairs and A, B, and C are 
real constants. The terms with B and C cannot cause scattering, since they are proportional 
t° Y2jS Sj ' Sj+S an d commute with the magnetic Hamiltonian. For the polarization picking 
up the A\ g mode, ej = (x + y)/\/2, ej = (x + y)/\/2, no intensity comes out by the same 
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reason (x and y are unit vectors pointing to the x and y axes). For the polarization picking 
up the Big mode, = (x + y)/y/2, e/ = (x — y)/y/2, we have 

H R .l^(S r S ; , x . S ; -S ; , y ). (A2) 

3 

Expanding this in terms of magnon operators, we obtain 

H R = ASzY^ (al(3^ + H.c. + ••■), (A3) 
k £k V J 

where 7^ = {cosk x — cosk y )/2. The scattering intensity I{oj) from this interaction is given 
by 

I(u) oc {ASzf V ( 2^ ) (-2)ImF(q = 0, u; k, k'). (A4) 

We calculate F(q = 0, u; k, k') by summing up the ladder diagrams including the Oguchi 
correction within the present formalism. We obtain the spectrum identical to that of Canali 
and Girvin, which is formed by a single peak at uj = 3.38 for S = 1/2.— 
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